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^ ! Abstract 

The Galilean invariance of the Navier-Stokes equation is shown to be akin to a global gauge 
symmetry familiar from quantum field theory. This symmetry leads to a multiple counting of 
^. infinitely many inertial reference frames in the path integral approach to randomly stirred fiuids. 

r~| , This problem is solved by fixing the gauge, i.e., singling out one reference frame. The gauge fixed 

theory has an underlying Becchi-Rouet-Stora (BRS) symmetry which leads to the Ward identity 
relating the exact inverse response and vertex functions. This identification of Galilean invariance 
\^ ', as a gauge symmetry is explored in detail, for different gauge choices and by performing a rigorous 

examination of a discretized version of the theory. The Navier-Stokes equation is also invariant 
under arbitrary rectilinear frame accelerations, known as extended Galilean invariance (EGI). We 
gauge fix this extended symmetry and derive the generalized Ward identity that follows from the 
BRS invariance of the gauge-fixed theory. This new Ward identity reduces to the standard one 
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' in the limit of zero acceleration. This gauge-fixing approach unambiguously shows that Galilean 

invariance and EGI constrain only the zero mode of the vertex but none of the higher wavenumber 
modes. 
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I. INTRODUCTION 



The formulation of the Navier Stokes equation with random forcing as a classical stochas- 
tic field theory l|, fl J Q, B, fl E] opens up the way to apply the methods originally developed 
for quantum fields Sj. This fact has been used in the majority of cases for applying renormal- 
ization group (RG) methods to models of fully developed turbulence This field theory 
approach has also been used on many occasions to derive exact relations between different 
correlation functions implied by Galilean invariance. These are akin to the Ward-Takahashi 
identities ( WTl) of quantum field theory lOl] . The most well known of these exact identities 
relates the vertex and response functions [2] • Based on this relation, claims have been made 
concerning the non-renormalization, under the renormalization group, of the advective or 
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13| . Despite the long history 



inertial term in the Navier-Stokes equation 
of applying these field theoretic methods, two major problems connected with the Galilean 
invariance of the path integral framework of randomly stirred fiuids have been brought to 
light recently. The first problem has to do with the functional itself. The standard dynamic 
functional for the randomly forced Navier-Stokes equation leads to spurious relations for the 
correlation functions involving the zero-mode of the fiuid velocity. This dynamic functional 
is therefore ill-defined. As shown in jl^ the dynamic functional for the randomly stirred 
Navier Stokes equation can be regarded as a gauge theory, and then the problem of the 
spurious relations is solved by fixing the gauge, which thus breaks the Galilean invariance. 
This is a new and important observation and leads to a well-defined functional expression, 
which has been missing for the past three decades. However, once the gauge is fixed, it 
is crucial to verify whether Galilean invariance (Gl) can be restored, and this is achieved 
by recognizing that the gauge- fixed theory is BRS invariant. This brings us to the second 
problem which has to do with the past mis-interpretations of the Ward identity that follows 
from the functional. This BRS invariance leads to a Ward identity whose physical conse- 
quence is that only the zero mode part of the full vertex, but none of the higher spatial 
modes, is constrained. In other words, contributions from the non-zero modes to the vertex 
can receive corrections under renormalization. This conclusion therefore provides greater 



clarity to, and in some cases goes against, previous interpretations regarding t 
consequences of Gl for the vertex in the Navier-Stokes equation 



le physical 



The present work supports the conclusions drawn in j?] in 12| and in 13|], and develops 
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in full detail the results presented in [13|, while also generalizing that work. For this, the 
gauge fixing concept for the NSE is explored in detail including the treatment of distinct 
gauge choices and the associated implications for the BRS symmetry left in the theory. 
Also, the analogy between the Ward identities found here for the NSE and those in QED is 
carefully examined. Finally, the considerations developed so far for Galilean invariance are 
thoroughly worked out here for extended Galilean invariance, which is a more general global 
symmetry of the NSE. We derive the consequences of this higher symmetry which leads to a 
new Ward identity for the forced NSE. The Ward identity for EGI has not to our knowledge 
been derived nor applied to the Navier-Stokes equation before. We then use it to see how 
the convective term in the Navier Stokes equation can renormalize. This has direct physical 
consequences. For this case, the time dependent but spatial zero mode part of the vertex is 
related to viscosity and mass renormalizations, while none of the higher spatial modes are 
constrained by EGI. These facts constitute the physical results of this paper. 

In order to consider the consequences of breaking Galilean invariance via gauge fixing, we 
start with the ensemble of histories of Navier-Stokes dynamics, and make direct use of the 



standard functional integral methods. The dynamic generating functiona' 



stirred incompressible fiuid has been extensive 



y studied for many years 
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for a random' 



[ii,y,i3i,y,B,y,j7 



16l . |17] to classical statistical dynamics [Ji 



It is based on the path integral approach [Ij, 
and is given by 

Z = J [DV][DcT]exp{-S[y,a]}, (1) 

where the action is ^ 

S[Y,(t] = dyidt J c/y(Ti(x,i(;)Aj(x - y)(Tj(y,t) 

-^ J d^dt + P..(V)^ - .oW,). (2) 

The instantaneous fiuid velocity is V, the conjugate field is cr; Dij{x) is the only non- 
vanishing cumulant of the stationary stirring force and x is a coordinate vector in d- 
dimensions. The projection operator is Pij(V) = {Sij — Vj^^Vj) and denotes the bare 



^ A Martin-Siggia-Rose type action can be derived for general Langevin equations of motion 19], and this 
includes the Navier-Stokes equation (NSE) with random forcing. The Jacobian that arises in passing from 
the general stochastic equation to thepath integral can, in effect, always be set to unity. This is treated 
in detail in Ref. 19 1. See also Ref. [l3| and Refs. 0,(1] which prove that this Jacobian is unity for the 
particular case of the forced NSE. 
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(unrenormalized) fluid viscosity. Now consider a second primed frame moving witli a con- 
stant velocity c witli respect to tfie unprimed frame. Then the relations between the in- 
stantaneous fluid velocity, conjugate field, time and coordinates of events in both frames are 
given by: 

V^(x,t) = V'(x',0 + c, (3) 
a%^,t) = ^'(x',t'), (4) 
t = t\ (5) 
X = X +ct. (6) 

"V^ denotes the result of the Galilean transformation applied to the velocity field V, and 
similarly for cr"^. Due to the Galilean invariance of the measure [DV*^] [-Dcr*^] = [-DV] [Dcr'] ^ 
and the action S under the transformation in Eqs.(j3}ED, is clear that the generating func- 
tional Eq.([T]) sums over all fluid velocity configurations. This results in multiple counting 
of physically equivalent configurations: namely, those that are equivalent up to a Galilean 
transformation (GT). In other words, this integral includes the sum over all inertial refer- 
ence frames. This feature leads to spurious relations among velocity correlation functions 
which must be removed in order to render a well defined generating functional for the theory 



13l |. The solution to this problem was demonstrated to be provided by gauge fixing [13 1. 
Moreover, if one wants to use the path integral for the Navier-Stokes equation for computing 
non-Galilean invariant quantities, like n-point velocity correlation functions, it is essential 
to carry out a gauge fixing along the lines we present here, and this is especially true for 



numerical implementations iJ]. Gauge fixing for Z can be viewed operationally in exactly 
the same way as in quantum field theory. In the present context, we must divide velocity 
configuration space into equivalence classes called the "orbits" of the GT. An orbit of the 
GT includes all velocity field configurations which result when all possible GT's are applied 
to a given initial velocity field configuration. Thus, Z is proportional to the infinite vol- 
ume of these orbits, and this volume factor is extracted out before defining this functional. 
However, once the Galilean invariance is explicitly broken, it is crucially important to verify 



^ The change of variables Eq.Q in the functional measure generates a unit jacobian: J = det ay^(y^'t'') — 

det ^^^gy^(^y, = det{6ij6{x' — y')6{T' — t')). Similarly for the transformation Eq.Q of the conjugate 
field. The Jacobian is also unity for the extended Galilean transformation in Eqs. (|44j|?7|) and Eq. lH^l) . 
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whether the exact identities relating inverse response and vertex functions continue to hold. 
We demonstrate that this is the case, because the resultant gauge fixed functional has an 
underlying Becchi-Rouet-Stora (BRS) symmetry. This symmetry leads to exact identities 
structurally similar to the Slavnov- Taylor (ST) identities of quantum field theory 2^. Here, 
the Slavnov- Taylor identity is significant because the ghost degrees of freedom make abso- 
lutely explicit that the vertex nonrenormalization only affects the zero mode. Contrary to 
previous assertions and claims, it is not even in the k limit that the claimed vertex 
nonrenormalization holds, but rather only at k = 0, the zero-mode. This is a subtle but 
very important point that must be born in mind. 

The outline of this paper is as follows. In Sec. Ill Al we pick a gauge which singles out a 
unique inertial reference frame, and carry out the gauge fixing procedure for the dynamic 
functional in Eq. ([1]). This leads to a gauge fixed action. In Sec. IllBI we write down the 
Becchi-Rouet-Stora (BRS) transformation that leaves this gauge fixed action invariant. This 
BRS invariance leads to the Slavnov- Taylor identities for the dynamic functional. These in 
turn imply the crucial Ward identity for the effective action that has been derived previously 



on numerous occasions, but always derived for the non-gauge fixed functional 



a 



a 



In Sec III CI we prove that this crucial identity remains valid, even after explicitly 



breaking the Galilean invariance, and selecting a fixed reference frame. Other gauge choices 
are possible, as illustrated in Sec III Dl The stochastic field theory for the NSE can be 
regularized in a space-time box and a corresponding Ward identity is derived in Sec III El 
This is used to demonstrate that the spurious relations first found in |l3'] for the continuum 
field theory arise in the regularized theory as well, and are therefore not an artifact of the 
continuum limit. The Navier-Stokes equation is also invariant under arbitrary rectilinear 
frame accelerations, a symmetry known as extended Galilean invariance (EGI) [2]]]. This 
is reviewed briefiy in Sec IIII Al Gauge fixing is carried out in Sec IIIIBI by choosing a 
single accelerating frame and the BRS transformation that leaves the corresponding gauge 
fixed action invariant is written down in Sec IIII CI To complete the analysis, we derive the 
Slavnov- Taylor identities for EGI in Sec IIIIDI and show that these lead to a new identity 
for the inverse response and vertex functions. This identity reduces to the standard one 
in the limit of vanishing reference frame accelerations. The Ward identities that follow 
from the EGI and GI in stochastic fiows bear a striking resemblance to the well known 
QED Ward identities that follow from the local U(l) gauge invariance. The similarities 
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as well as the differences are spelled out in Sec HVl We discuss the implications of both 
Galilean and extended Galilean invariance for vertex renormalization in Sec |V] and present 
our conclusions in Sec IVII The main steps for gauge fixing the Galilean invariant dynamic 
functional are collected in an Appendix which closely parallel the standard technique as 
employed in quantum field theory. 



II. GALILEAN INVARIANCE: GAUGE FIXING AND BECCHI-ROUET-STORA 
SYMMETRY 

A. The Vo^ gauge 

For the reasons stated above, we proceed to gauge- fix the dynamic functional in Eq.([T]). 
The basic outline for doing so in an arbitrary gauge f(V) is summarized in general terms 
in the Appendix [Al This is intended to parallel closely the procedure as it is used in 
quantum gauge theories. Once we choose the specific gauge introduced below, we then only 
have to calculate the corresponding Fadeev-Popov determinant Eq. flASp . and then select the 
arbitrary function U that appears in Eq. (1A6l) to complete the procedure. Here and below, 
we consider spatially bounded fluid systems. 

Recall that the purpose of gauge fixing here is to single out one inertial reference frame. 
This can be accomplished by constraining the zero mode part of the full instantaneous 
velocity field, because this mode corresponds to the constant velocity of the system as a 
whole, that is, its bulk velocity. We thus make the gauge choice specified by 

f(V) = ! cixdtV(x,t), (7) 



vol _ 

where vol = L^T is the (finite) volume of the space-time box bounding the fluid. This 
projects out the zero mode from the full instantaneous velocity field V, as required. 
Next, from Eqs. (l3|5l[6l) . we have 

V^^(x,t) = V,'(x-ct,t) + c,-, 

dV'ix, t) 

= V^'(x, t) - Ckt + cj + 0{c'), 



5v^^(x,t) = v;-vj' = ckt 



Thus ' ^ Lc=o = t ' ^ - 5^^. (8) 



Then from Eqs. (17f8|) we have 



— ^ 'c=o = -7 / d^dt (t^^ - 6,, ) . (9) 



dcjJ vol J \ dxj 

This simphfies due to boundary conditions. Imposing the physically reasonable periodic 
boundary conditions for a spatial box of side length L, we have 

dV(:si t) 

J "^^i " = ^' ^) ~ = 0' ^) = 0' (10) 

where are the coordinates orthogonal to Xj. Thus for the Fadeev- Popov determinant in 
Eq. llASp . we calculate 

A/[V] = det (1^) |e=o = det(-5,,) = j dridrj* exp - zry* • r/). (11) 

The last equality expresses this determinant as an integral over constant complex conjugate 
Grassmann vectors rj and rj*, with {rii,ri*} = fiol]. Lastly, since our gauge f(V) Eq.d?]) 
is simply a constant vector, we repeat the arguments leading to Eq. (]A6P and Eq. (lASP in 
the Appendix, but here replace the function U(x, t) ^ b by a constant vector in Eq. flA7p . 
and drop the spacetime integrations indicated there, that is, instead of Eq. (1A7P we choose 
G[h] = exp ^ — j , where ^ > is a free parameter. 

Putting these results together, and after integrating over the delta -function in Eq. flA6p . 
which in this case is S(Vo — b), the gauge fixed action reads: 

SGF = S + ^Vo' + ^V*■r^, (12) 

where S is given in Eq.(l2]) and we define Vq = ^ / c?X(it V(x, t), and note that this has the 
correct dimensions of a velocity. Galilean invariance is manifestly broken in Sgf- 



B. Becchi-Rouet-Stora invariance of the gauge-fixed action 

Despite the fact that Galilean invariance has been broken explicitly, the gauge-fixed 
theory Eq. fll2l) does have a fundamental underlying symmetry which makes full use of the 
Grassmann vectors introduced to represent the Fadeev- Popov determinant Eq. flTT]) . This 
symmetry was originally established to be a general consequence of gauge fixing in non- 
Abelian gauge theories, and is known as Becchi-Rouet-Stora (BRS) invariance^ j20|. 

^ The BRS symmetry in our paper arises from gauge-fixing, and should thus be distinguished clearly from 
the BRS symmetry that arises from regarding the stochastic differential equation itself as the constraint; 
see, e.g., Ref. 0|. 
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Consider the following infinitesimal BRS transformation where (^^ = is a real Grassmann 
constant: 



5brsx = -({ri + vl, (13) 
SBRst = 0, (14) 
5BRsV(x,t) = CiVk + vDi^^^^-Ciri + ri*), (15) 

<5brs^(x, t) = dm + Vl)t^^^, (16) 

SbrsV = -^VoC, (17) 

^BRsr;* = +^VoC. (18) 

Note that Eg. (1151) automatically implies the corresponding BRS transformation for the ve- 
locity zero mode Vq: 

WVo = ^ y c?xc/t WV(x,t) = -C(r7 + rj*). (19) 

Now Eqs.( fT3lfT6l) are recognized to be just the infinitesimal version of the Galilean trans- 
formation Eqs.(l3H6]), but with the boost velocity replaced by the Grassmann vectors 
c —>■ + ^*)- Therefore, the action S and path integral measure are automatically 
invariant under this subset of the full BRS transformation. Next, from Eqs. fll7lll8|[T9!) . we 
easily verify that the combined gauge fixing and Grassmann terms in Sgf are invariant. 
From these observations we therefore conclude that the complete gauge fixed action Sgf in 
Eq.( fT2l) is invariant under the full BRS transformation Eqs. (11311181) : 

SbrsSgf = 0. (20) 
This BRS invariance leads to exact identities which we derive below. 



C. The Slavnov- Taylor and Ward identities 

We introduce source terms into the gauge-fixed functional as follows: 

Zgf[3, S, d, e*] = J [DV][D(T]dr]dr]* exp{-^GF[V, cT,r],rj*] + 6* ■ + ■ t]* 

+ y"rfxc/t{J(x,t) ■ V(x,t) + S(x,t) ■o-(x,t)}}, (21) 



where d and 6* are complex Grassmann vectors. Next, we subject this functional to the com- 
plete BRS transformation (that is, we displace all fields by the infinitesimal BRS transfor- 
mation) Eqs. fll3tiI5]) . Since the measure and the gauge- fixed action Sgf are BRS-invariant, 
only the source terms in Eq. fl^T]) will be affected. Moreover, as = 0, we can easily expand 
the exponential: exp((^yl) = 1 + C,A. The gauge-fixed functional Eq. fl21l) thus transforms as 
Zgf — ^ Zgf + ^brsZgf, where SbrsZgf = can be written as follows: 



C vol J 5J(x,t) ^dOj oey \ 



ZGF[J,S,0,r] =0. (22) 



Here the operator 



(5 b 

Oj = I d^dt ( JmtVj—- + T.mtVj-— - Jj ) (23) 

J rn Z^rn. 



is the same operator appearing in the well known Ward^dentity that follows from the 

^: that is, 6jZ[J, S] = 0. 



Galilean invariance of the non-gauge fixed functional 

To finish, we introduce the effective action F and write the Slavnov- Taylor identity Eq. fl22|) 
in terms of this quantity. Let W = In Zgf, then the generating functional of one-particle 
irreducible functions is given by the Legendre transform 

r[V,z, a,i, ri,i, ri2 = -W[3, S, 9, 6*] + 0* . ry,, + 6> • ry* + J dkduj (J ■ V,, + S ■ cr,;), (24) 

which for convenience, we express in wavevector k and frequency space u. Here, the label 
"c/" reminds us that the fields so indicated are averaged over the fiuctuating fields and in 
the presence of the source terms: that is, using the ensemble in Eq. (l2T|) . From Eq. (!2^ we 
have that 

We can now straightforwardly write down the identity for F that follows directly from 
^brsZgf = 0: 

/ar dT 



du 5V^{Ku) ' du 5(T^(k,cj) ' ' ' '6Vf(k,u) 

We now make use of this formula and apply it to the problem at hand. The dependence of 
F on 77 and 77* is simple since these constant Grassmann vector fields do not interact nor do 
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they couple to the velocity or conjugate fields ^. Thus, we can immediately write 



r [Yd , (Tci , r]ci , ri*i] = -iri*i -rid + T [Y^ ,(t,i], (27) 

where r[Vci,o"c/] does not depend on either rjd or r]*i, and the first few required terms are 
written out in Eq. flBip in Appendix [Bl 

The pertinent identity we seek is then obtained by inserting Eq. fl27l) into Eq. (l26|) . and 
then differentiating this with respect to 6/6Vf''(k,uj)6/6a^{—]<.,—uj), and then setting V^^ 
= cr'^' = 0. Note that the terms in Eq. fl26p depending on the gauge parameter ^ do not 
contribute to this sequence of steps, and most importantly, we end up obtaining the following 
result: 

iVci + v:i)j {''^g^^in'^ (k, -k, -c^) + rSt^ (0, 0; k, u- -k, -u;)) = 0. (28) 

Moreover, since rjd and r]*i are arbitrary, the expression within the larger parentheses must 
vanish identically. T^^'^^ and F*^^'^^ denote the inverse response and vertex functions, respec- 

cnown Ward identity derived previously on many occasions 



tively. Thus, we recover the well- 
from the non-gauge fixed action 

mm 

. Alternatively, we can differentiate Eq. (!26|) with 
respect to d/d7]f 5/5 V^'^' (k, tu) 5/(5(7^' (—k, —uj) followed by setting 77*^' = 77*; = V^' = cr*^' = 0. 
The end result is the same. 

In summary, the Ward identity Eq.( l28i) relating the exact inverse response function to 
the exact vertex of the NSE holds for the gauge fixed action Eq.( JT2l) . is a direct consequence 
of the BRS invariance of the latter, and is independent of the gauge parameter ^. From 
this point on, the entire discussion regarding the implications of Galilean invariance for the 
vertex renormalization [tI] continues to apply; see Sec |Vl Most importantly, we see that 
those arguments are valid even after breaking the Galilean invariance and going to a fixed 
reference frame. 



D. The (5- gauge 

Other gauge choices are possible. To briefiy illustrate this, we consider this field theory 
in the 5-gauge, in which we work directly with the functional integral in the form Eq. (1A6I1 



^ This dependence of V on the Grassmann variables is similar to that for gauge-fixed quantum electrody- 
namics, see, e.g. the development in Ramond's book 
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(and with the overall infinite term removed). Put f(V) = Vq and U = b in Eq. (lA6|) . choose 
as the weight function G[\J] = 1 rather than Eq. flATp . and use the representation 

S'iVo -h) = j exp [^K . (Vo - b)] (29) 



^BRsVoj 






= K,C, 


SbrsVj 


= -K,C, 


^BRsb 


= (5brsK = 



In this case, by using the final equality in Eq. ffTT]) . for expressing the determinant of the 
identity matrix ClS db Grassmann integral, it gives 

Sgf = S[V, a]-iK- (Vo - b) + irj* ■ rj. (30) 

The Galilean transformations are as above Eqs.(l3]- [6l), for which ^[V, cr] and the measure 
are invariant. For the BRS transformation, we identify c Ci'H* + ^) before with 

(31) 
(32) 
(33) 
(34) 

We point out that 77, r}* need not be a complex conjugate pair; see footnote in Sec IIII Ci 
With this, one can check 

Sbrs {^V* . r7 - zK ■ (Vo - b)) = 0, (35) 
proving that the gauge-fixed action in this gauge Eq. (l30l) is BRS invariant: SbrsSgf = 0. 

E. Ward identity in a box and spurious relations 

^n order to remove any ambiguity from our treatment of the continuum NSE path integral 
in 3] and the spurious relations we found there, we here perform a similar analysis but 
now for a discrete and hence regularized version of NSE field theory. In this Section we 
use the dynamic functional Z for the NSE derived in a box (see Appendix [C]) to obtain the 
associated Ward identities. An important technical distinction stems from the fact that in a 
box, there is a smallest limiting finite velocity boost. Only in the continuum can infinitesimal 
boosts be treated. 
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To begin, and including the sources from the outset, we start with 
Z[J,S] = / Yldva{n,j)daa{n,j) X 

exp - S[v,cr] + ^ (Jai-m, -l)vaini,l) + S^(-m, -/)(Ja(m, (36) 



Here, the discrete action S is given by Eq. flCQp . We repeat all the steps that for the con- 
tinuum case, lead one to obtain the Ward as for example in 7|]. Here, we carry this out 
for the discrete theory. So we next transform Z using the discrete version of the Galilean 
transformation Eq. flC12p . Doing so we obtain: 



J Yl d-Vain,j)daain,j) exp ( - 5'[v, cr] + y] ( Ja(-m, -l)va{m, I + c-m) 

a,n,jr m,Z 

c ■ J(-m, -l)6^{in)6{l) + S«(-m, -Z)(Ta(m, j + c-m))). (37) 



Now unlike the continuum case, we cannot consider an infinitesimal boost velocity. Instead, 
the best we can do is consider the smallest nonzero boost velocity, which, according to 
comments immediately below Eq.( lC12l) . is 6c = (1,1,1). We could as well have chosen 
(1,0,0), or (0,1,0), or (0,0,1), or (1,1,0), etc. The point is, the smallest allowable boost 
velocity is not infinitesimal. 

We next emulate what is done for the continuum case, namely Taylor expand the velocity 

n 

and conjugate fields about zero boost, see for example, Eq.(4) in [7]. The best we can do is 
to employ the following approximation: 

^^(m, / + c ■ m) = t;Q,(m, /) + c ■ mAt;Q,(m, /) 

= VainiJ) + {mi + 1712 + ■ms)AvaimJ), (38) 

where the finite difference operator A/(m, /) = /(m, / + 1) — /(m, /). 

To finish the derivation, we insert these finite difference approximations back into Eq. (!37Il . 
This implies that the functional transforms as Z ^ Z + 6Z. We require that the extra terms 
proportional to the boost velocity c vanish identically. This leads to: 

d 



Ja(-m, -/)mjA 



dJa{—m, —I) 



+ S«(-m, -l)miA-^^j^^^^ - J,(-m, -/)5^,o<^«,o )Z[3, S] = 0. (39) 
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This is the Ward identity for Z and should be compared to the Fourier transform of the 
continuum version Eq.(6) in j?]. 

Differentiate this WTI for Z, Eg. (15^ . with respect to aj(„m -/) I J=s=0) to obtain 



mi A < fA;(m, /) > —b, 



ik 



0. 



(40) 



This is well behaved and never singular, but it suffers from the same defect as its continuum 
version [is!] when we set i = k and sum over this index: fluid incompressibility implies 
mifi(m, Z) = 0, which since it is zero, leads to a contradiction in Eq. lHOl) . 

For higher point correlations, more spurious relations are found from Eg. (139!) . similar to 



13|. The source of these spurious relations is the same 



those found in the continuum case in 

as in [l3|] , and we can verify that the functional integral in Eq. (1371) has an infinite prefactor 
due to the integration over the velocity zero mode. This can be fixed by the same gauge 
fixing procedure as we carried out in ISj, taking appropriate modifications for the discrete 
case. Most importantly, the treatment in this Section demonstrates that these spurious 
relations are not an artifact of the continuum path integral expressions. 



III. EXTENDED GALILEAN INVARIANCE (EGI): GAUGE FIXING AND BRS 
SYMMETRY 



Galilean invariance admits an interesting physical extension. As it turns out, the Navier- 
Stokes (NSE) equation is invariant under rectilinear frame accelerations, a symmetry known 
as extended Galilean invariance (EGI) 211]. This invariance was noted and used in the 
context of the KPZ equation 221], as well as in the stochastic Burgers equation [23], it has 



24j . and finds practical 



been used in probability density functions in turbulence modeling 

applications in fluid animation simulations [25^. We review this invariance of the NSE 
below. We then proceed to gauge-fix this symmetry, write down the corresponding BRS 
transformation that leaves the gauge- fixed action invariant, and then deduce the associated 
Slavnov- Taylor identity for the functional and the new generalized Ward identity that follows 
from it. 
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A. Extended Galilean invariance of the Navier-Stokes equation 



Consider the Navier-Stokes equation (NSE) in an inertial frame E: 

Ot OXj OXi 

where u is the kinematic viscosity of the fluid, Vj(x, t) and n(x, t) are the instantaneous 
values of the velocity and pressure, and the continuity equation takes the form 

for an incompressible fluid. In this case, the density is constant for all x and t so for 
convenience we may work in a system of units where the fluid density is taken to be unity. 
Also, as is well known, taking the divergence of each term in Eq.l HTl) . and invoking Eq.( H2il . 
leads to a Poisson-type equation for the pressure, viz., 

v^n ^ (43) 

and this result is useful for establishing the Galilean transformation of the fluid pressure. 

Now consider a second noninertial frame E moving with respect to E with a variable 
but rectilinear velocity c(t). The transformation between these two coordinate systems is 
as follows, 

x' = x-A(t), (44) 
t' = t, (45) 
VV,0 = V(x,t)-A(t), (46) 
n'(x',t') = n(x,t)+x'- A(t), (47) 

where A(t) = c{s)ds, and the dots stand for the derivative taken with respect to time. 
Applying these transformation rules to (HTl) proves that the transformed NSE takes the form 

which is the NSE in the noninertial accelerating frame. Note that X{t) = c{t) = A(t), so 
that the frame acceleration is absorbed into a modified pressure Eq.f l47l) . This demonstrates 
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that the NSE Eg. (1411) is invariant under rectihnear frame accelerations. It is clear that when 
the above transformation rules Eqs. (I44ti^ include one for the conjugate field: 

cT'i^',t')=ai^,t), (49) 

then the action Eq.([2]) and dynamic functional Eq.([T]) are invariant as well under this ex- 
tended Galilean transformation. 



B. Gauge Fixing for Extended Galilean invariance 

We next gauge- fix the EGl of the dynamic functional in Eq.([T]). Once again, we refer to 
the general gauge-fixing procedure as outlined in Appendix |Al We make the specific gauge 
choice introduced below, calculate the corresponding Fadeev-Popov determinant Eq. (lA5p . 
and once again select the arbitrary function U that appears in Eq. flA6p to complete the 
procedure. 

In the case of EGI, the purpose of gauge fixing is to single out one rectilinearly accelerating 
reference frame. This can be accomplished by now constraining the time-dependent zero 
mode of the full instantaneous velocity field. This mode corresponds to the acceleration of 
the system as a whole, that is, its bulk acceleration. We therefore make the following choice 
for the gauge function, where vol = is the volume of the spatially bounded box: 

f(V) = ^ y f/xV(x,t) = Vo(t). (50) 

This projects out the instantaneous and arbitrary time dependent bulk velocity of the entire 
bounded system as required. From Eqs. (I44tii6|) we have 

\A^x,t) = V^'(x',t') + Ai(t), (51) 
= V^{^-X{t),t) + X,{t), (52) 
= y'(x,t)-Afc(t)^V;'(x,t) + A,(t) + 0(A2), (53) 



k 



dx 

8\m{t') ~° V dXm ^"^dt^ 



-^6it-t'). (54) 

The FP determinant Eq. ( 1A5I) in this case works out to be, again using periodic boundary 
conditions for the bounded spatial box. 



J Vr]Vr]*exp [- ^ J dtrjit)* ■ j^nit))- (55) 
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The functional integral in Eg. (1551) is over two time-dependent Grassman vector fields r){t) 
and ?7*(t). For the convergent integral over the gauge group volume, we now choose (see 
Eq.dMl) in Appendix E]) 

G[b] = exp(-iy" dth^t)^ (56) 
which yields the gauge-fixed action for EGI: 

Sgf = S[V,a] + j dt |lvo^(t) +zr7*(t)^r/(t)| . (57) 

At this point, Extended Galilean invariance is manifestly broken in Sgf- 



C. Becchi-Rouet-Stora transformation for gauge fixed EGI 

We begin by writing down the BRS transformation. We observe that the transforma- 
tions in Eqs. (I58tl62|) below are just the infinitesimal version of the extended Galilean trans- 
formations Eqs. (l44ti47|) . but with the variable rectilinear boost velocity now replaced by 
X{t) = c{t) —* (rjit), where C is a real constant Grassmann parameter: ("^ = 0. Therefore 
the nongauge-fixed action S and the path integral measure are automatically invariant under 
this subset of BRS transformations: 

(58) 
(59) 
(60) 

(61) 

(62) 



^BRSX 


= -Crjit), 




^BRst 


= 0, 




5BRsV(x,t) 


. ,.^^V(x,t) 


-cm, 


'5BRSO-(x,t) 


(9cr(x,t) 




5BRsn(x,t) 


^ C..(t)'^i^'^^ 


+ X • cm 



To complete the BRS transformation, we need to know how the velocity zero mode Vo(t) 
and Grassmann sector transform. The remaining BRS transformations are given as follows: 

5BRsVo(t) = -cm, (63) 
SBRsm = 0, (64) 
SBRsr]%t) = yVo(t)C. (65) 
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Note that Eq. (l63i) is an immediate consequence of Eq. (160!) . Secondly, note the rather different 
ways that r] and 77* transform. This "asymmetry" is in fact famihar from non-abehan 
quantum field theory ^ 

Summarizing up to this point, we have demonstrated that the gauge fixed action Eq. flFT]) 
is invariant under the above BRS transformation: 

SbrsSgf = 0. (66) 



D. The Generalized Slavnov- Taylor and Ward identities 

To conclude, we subject Zgf, which we now define below, to the BRS transformation 
and derive the generalized Ward identity that holds for F, 

ZGF[J,^,0,e*] = j[D\][Da][Dri][Dri*]exp{-SGF[y,(T,ri,ri*] + j dt {6* ■ + ■ rj*) 

+ ycixdt{J(x,t)- V(x,t) + S(x,t)-<T(x,t)}}. (67) 

Here of course, the gauge fixed action Sqf is given by Eq. (l57|) . and the Grassmann vector 
sources 0{t), 0*{t) are now time dependent. We next subject this functional to the complete 
BRS transformation listed in Eqs. (l58H65l) . Just as in the case treated in Sec III C\ the 
measure and gauge-fixed action are invariant and only the source terms will be affected. 
The remainder of the steps involved are similar to those employed above. We expand the 
exponential to first order in (, and then make use of the Legendre transform similar to 
Eq. (1241) to define the effective action F. For the present case, the Grassmann source terms 
are integrated over time in the Legendre transform. After a few simple steps we arrive at 
the following Slavnov- Taylor identity: 



0. 

(68) 



^ See for example the chapter on "Yang-MiUs Theory: Slavnov- Taylor Identities" , page 370 in [l^ for full 
details of the proof of the BRS invariance of gauge-fixing and ghost terms in a general Yang-Mills gauge 
theory. Note also the use there of two independent real Grassmann fields. The limit of vanishing Lie 
algebra structure constants (the Abelian limit) in Yang-Mills quantum field theory implies that one of 
the two ghost fields transforms to zero under BRS. See also remarks [3] to the effect that rj and tj* can 
be independent, it is not necessary to think of them as hermitian conjugates, indeed they can be any two 
independent anticommuting fields. 
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The pertinent Ward iderrtity we seek is obtairred by inserting the functional Taylor series 

for r (see e.g., Eq.(9,10) of Ref.pl]) into Eq. fl68l) . differentiating the latter with respect to 
Vfc(y,t') and cr(w,t"), and then setting all the fields equal to zero. We immediately see 
that term depending on the gauge parameter ^ will not survive this sequence of steps and 
furthermore, since 77^' (t) is an arbitrary function of time, the identity Eq. fl68l) effectively 
reduces to 

/rf^x fVjV;(x,t) , + Vj(Tf(x,t) ^ ^ + e-x//^ =0. (69) 

Carrying out the specific differentiations mentioned above yields a relation between the exact 
inverse response F*^^'^^ and the vertex functions F*^^'^); 

I / rf'^xFSj)(x,t;y,t';w,t") = 6{t - t")^jif {y,t';w,t) 

+ 6{t~t')^T^f{y,t;wX). (70) 

This is the new Ward identity that follows from EGI in configuration space and time, derived 
here for the first time. Whereas the vertex function in the conventional Ward identity (see, 

n 

e.g., Eq.(ll) of Ref.|'il]) is integrated over both coordinates and time, here by contrast, 
the vertex is integrated only over coordinates and is differentiated with respect to time. 
Using translational invariance (in space and in time), Fourier transforming ( 170|) and then 
integrating over J dt J dt' J dt" e~*'^i*e~*'^^*'e~*'^^*" to get rid of the delta-functions, we obtain 
the Ward identity in wavevector and frequency space (note: here, u stands for a frequency, 
not the viscosity): 

knT^kj^\k,(j + iy-,-k,-uj-iy)-knT'j^.^\k,iy-,-k,-iy) = -ujT^^^f {0,uj;k,iy; -k,uj -u). (71) 
If we take the zero frequency limit 0, we then obtain the identity 

kn^rif{k,u--k,-u) = -FSj)(O,0;k,z.;-k,-z.), (72) 



which is precisely the Ward identity Eq. fl28|) that follows from the standard Galilean invari- 
ance of the NSE, that is, when the frame E moves with a constant velocity X{t) = with 
respect to the lab frame E [7|. 

For EGI, there is second way to implement the BRS transformation, intimately related 
to how we choose to define the gauge-function. Thus, if we now choose X(t) instead of X(t) 
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as the gauge function, then in place of Eq. (lMI) . we would have, after using the identity 

^^kit) = Jq du '^^'^^^ = du Xk{u), the relation 



-5,m5{t-t'). (73) 



The FP determinant, invoking the same boundary conditions as before, is now given by 



A^[V] = det(^^^)U=o = det(-5,^5(t-t')) 



= j VriVri*exp(^-i J dtr}{ty -riit)^, (74) 
and then finally Eg. (156!) yields the gauge fixed action 

Sgf = S[Y, ct] + Jdt {^Vo'(t) + iri*{t) ■ ri{t)^ , (75) 



to be contrasted with Eg. (1571) . The net result of this alternative choice of gauge function is 
to remove the time derivative from the Grassmann fields in the gauge fixed action Eg. (175!) . 

In this latter case, the BRS transformation that leaves the gauge- fixed action Eg.( 175|l 
invariant includes in part the infinitesimal extended Galilean transformation obtained by 
replacing the variable boost velocity by the time dependent Grassmann vector c{t) C'nit)^ 
where ( is the real Grassmann constant we introduced before. Namely, 



^BRsx = -C / dur]{u), (76) 
Jo 

SBRst = 0, (77) 

5BRsV(x,t) = ( lyurj.iu)^^^^ - Cr7(t), (78) 

(5BRso-(x,t) = dumiu)^^^^, (79) 

5BRsn(x,t) = C jyurikiu f^^^^K ^- (80) 

This version of the BRS transformation reduces to that found in 13| in the limit that EGI 
goes over to standard GI (zero frame acceleration). To complete the BRS transformation, 
we must specify how the velocity zero mode and Grassmann variables transform. These are 
simply 

5BRsVo(t) = -Crj{t), (81) 

SBRsriit) = 0, (82) 

SBRsri%t) = yVo(t)C. (83) 
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Note: we could also choose the more "symmetric" assignment c(t) C(^(^) + 'n{t))y ^^'i 
then make the associated simple adjustments to the above set of transformation rules, but 
this has no effect upon the generalized Ward identity that follows from BRS. To summarize 
so far, the gauge fixed action Eq. fl75l) is invariant under this BRS transformation Eqs fl76ti8"^ . 

To finish, we next subject Zgf in Eg. fl67p to the above BRS transformation, where now 
Sgf is given by Eq fl75l) . Following the now rather familiar steps, this leads to the equation 
^brsZgf = 0, and expressing this in terms of the effective action F, yields the Slavnov- Taylor 
identity 

. / 5F „ 5F „ d 6T 



dt 5Vf 



= 0, 
(84) 

where we have defined p'^\t) = du rj^^u), so that p^^{t) = ■q'^\t). The Slavnov- Taylor iden- 



tities that follow from EGI Eqs. (]68ll84p and the Ward identities that they imply Eqs. (l69ll7Tll 
have not, to our knowledge, been derived nor considered before. From this point on, the 
entire discussion immediately following Eq. (!68l) remains intact and applies without modifi- 
cation to Eq. (1511) , as does the discussion concerning the implications of this identity Eq. ([Hi 
for vertex renormalization. 



IV. WARD IDENTITIES IN QED AND IN THE STOCHASTIC FIELD THEORY 
OF RANDOMLY STIRRED FLUIDS 

The well-known Ward identity in QED relating the exact photon-electron vertex and the 
electron propagator reads, 

S'f'\p + q)- S'p-\p) = g%(p + q,p), (85) 

where p+q, q denote the four-momentum of the entering and emerging electron lines, respec- 
tively, F^ denotes the exact photon-fermion-fermion vertex and S'p{k) is the exact Feynman 
fermion propagator [27]. This identity allows us to compute 5*^ directly from knowledge of 
the vertex and is represented graphically in Eq. (l86l) where the blobs denote the exact ex- 
pressions. The four-momenta carried by the individual fermion lines (single arrowed-lines) 
is indicated by the quantities within parentheses. Note that the four-momenta carried by 
the photon line (wiggly) on the right-hand side of Eq. (l86|) is therefore = q. This identity is 
obviously satisfied by the bare vertex F^ = 7^ and bare propagator, Sf{p) = (p^7^ — m)^^ 
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(p+q) 




(86) 

Taking the g — limit of Eq. (185|) yields 

'-^-M (87) 

which holds for the special case when the momentum transfer to the photon approaches 
zero. The crux of this identity is the observation that the derivative of a fermion line with 
respect to the four-momentum is equivalent to the insertion of a photon line in the limit of 
zero momentum transfer. This derivative rule therefore generates the "soft photon" vertex 
in QED. We illustrate this differential identity Eq. fl87p diagramatically in Eq. fl88l) : note that 
the four-momenta of the photon line (wiggly) is zero: g = 0. 



^{ (P)} = (P) (P) (88) 




Having briefly reviewed the main features of the key Ward identity in QED [27|, we come 
back to the analogous identities that hold in the stochastic field theory of the randomly 
forced NSE. The exact identity Eq.l lTTI) is shown in diagram form in Eq.( l89l) and holds for 
the exact inverse response function and triple vertex (indicated with blobs). The zig-zag and 
arrowed double lines denote the conjugate and velocity fields, respectively. The momentum 
and frequency carried by the various lines is indicated by the associated pair of variables 
(fc, uj). Note that the momentum and frequency carried by the vertical velocity line (arrowed 



double line) on the right-hand side of Eq. fl89l) is therefore = (0, uj). This identity is obviously 



satisfied by the bare inverse response function T[^-'^\k,uj) = {—iuj + uok'^)Pij{k) and bare 



(2, 1) 

vertex Fj^j'- = iknPij{k) This identity is structurally quite similar to the QED Ward 
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identity in (l85l) but there is an important difference: in Eq. (17T]) the 3- momentum transfer to 
the inserted velocity hne is always zero, whereas any finite frequency u can be transferred. 




Taking the zero-frequency limit {uj 0) yields the further differential identity in Eq. (ff2ll . 
which as we see, holds for the special case when the frequency transfer to the inserted zero- 
momentum velocity line goes to zero. The crux of this identity is the simple observation 
that the derivative of the exact response function with respect to frequency is equivalent to 
the insertion of a velocity line in the limit of both vanishing momentum and frequency trans- 
fer: = (0,0). We illustrate this differential identity in Eq. (!90!) . This derivative operation 
therefore generates the "so/^ velocity" vertex in NSE diagrammatic perturbation theory: the 
momentum-frequency pair carried by the velocity (the vertical arrowed double line on the 
right-hand side of Eq.(l90])) is = (0, 0). 




The rigid frame velocity, relative to the lab frame, is what enters into the WTI's for the 
NSE. This velocity can either be a constant, in which case we have the strict Galilean trans- 
formation, or it can be time dependent but rectilinear, in which case we have the extended 
Galilean transformation. This velocity enters the triple vertex as a background field. This 



is further supported by looking directly at the NSE in momentum space and time 

22 



28| and 



examining the nonlinear contribution. This contribution sums over all velocity modes and 
time dependences; there is a single contribution coming from the time dependent zero-mode 
which is simply the contribution coming from rigidly translating the reference frame. So the 
WTI is the statement that it is only this rigid motion that does not renormalize. This bulk 
system motion can correspond to a constant velocity, or a time dependent but rectilinear 
velocity, i.e., an arbitrary rectilinear acceleration. 

Regarding our comparison to gauge field theory, this has a two-fold motivation. On the 
one hand, we note that the fluid velocity transformation V ^ V = V + c and the transfor- 
mation of the QED vector potential A ^ A' = A + V0 are affine transformations leaving the 
convective derivative in the Navier-Stokes equation invariant and the U{1) covariant deriva- 
tive invariant, respectively. The second motivation is inspired in large part by the structure 
of the Ward identities that follow from breaking Galilean (or extended Galilean) invariance, 
since these identities Eqs.f l71|l72!) do closely resemble those that follow from breaking gauge 
invariance in QED, Eqs.( !85l[87l) . respectively, as we argued above. QED of course has a local 
(gauge) invariance, whereas Galilean invariance is a global symmetry. 

One can also make an analogy between Galilean invariance and field theories with global 
symmetries. For example consider symmetry breaking in a globally invariant field theory 
such as N-component scalar field theory with a global 0{N) invariance, where many results 
have been obtained. Considering this case in the symmetry broken phase and shifting the 
field about the vacuum expectation value v = (0), one can derive Ward identities involving 
V. In particular, there is an identity relating the two and three point IPI functions (vertices) 
that can be compared to the one we derived for broken GI or for broken EGI. This is (for 
complete details refer to IJ]), 

tu^fJip) + n^r'i;\p) + t;,t^,rS(p)(o,p, -p), (91) 

where the t", a = 1,2, ...,N are N x N real antisymmetric matrices: the generators of the 
0{N) Lie algebra. These of course have no analogue in our field theory, as the Galilean 
transformation is abelian, whereas 0{N) is a nonabelian symmetry group. Apart from 
these important Lie algebra factors, this identity holds for the triple vertex r*^^-* at zero four- 
momentum transfer, an important distinction to Eq.( 17Tl) which holds for zero 3-momentum 
transfer but for arbitrary finite frequency transfer. From the asymmetry of the we can 
write = — t"^ and so express the first two terms in Eq. flQTj) as a difference, but there is no 
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obvious limit to take here that would yield even a formal derivative, thus making comparison 
to the zero-momentum and zero frequency transfer Ward identity in Eq. (l72|) difficult. From 
these considerations, we conclude that the Ward identities resulting from broken Galilean 
and broken extended Galilean invariance bear a closer resemblance to those arising in QED, 
than to the those from the above globally invariant N-component scalar field theory. 



V. IMPLICATIONS FOR VERTEX RENORMALIZATION 

The vertex function F*^^'^^ is associated with the nonlinear or advective term in the Navier- 
Stokes equation. In the action formalism Eq.([2]), this vertex is given at tree- level by the 
nonlinear or advective term multiplied by the conjugate field: (TkPkji'V)^^^^- . It is thus a 
three-legged object, built up from one conjugate field and two velocity fields, and this fact 
is refiected when writing out the three arguments of wavevector-frequency pairs. Of course, 
both wavevector and frequency are independently conserved at the vertex, so it is in general 
a function of two independent wavevector-frequency pairs (see e.g., Ref.jJ] for the elements 
of Navier-Stokes diagrammatic perturbation theory). The most general form of the vertex 
taking into account corrections and possible renormalization effects must reflect this fact. 
We can separate out the tree level or zero loop contribution, and thus the term A represents 
all the possible (higher-loop) corrections: 

(2 1) 

Tlj'^'{q,Lj; k,u;q-k,uj - u) = ikiPjkik) + Aijk{q,uj; k, z/; q- k,u; - u). (92) 

In Eq. ( 192|) our convention is that the middle argument corresponds to the incoming wavevec- 
tor and frequency carried by the conjugate fleld, whereas the flrst and third arguments cor- 
respond to the wavevector /frequency pair carried by the two velocity flelds that meet at the 
vertex. 

In a similar fashion, the exact inverse response function can be written as follows: 

Tff:^\k,uj; -k, -uj) = [-iuj + u{uj, k)e + S(cu, k)] P,-fc(k), (93) 

where ^{oj, k) and S(u;, k) denote a renormalized viscosity and a "mass" term, respectively. 
At tree level, this corresponds to the term in the action cr^ — v^V'^Vk^ . So, the inverse 
response function F'^^'^-', is a two-point object (built from one coniugate and one velocity 
field) and thus a function of just one wavevector and one frequency ^. The Ward identities 
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derived above in Eqs. fl71ll72p impose certain constraints on these three functions z/, E and A. 
In other words, they relate the nonhnear and hnear parts of the Navier-Stokes equation. First 
consider imphcations of extended Gahlean invariance. Inserting Eqs. (l92f93p into Eq. (17T|) 
implies that 

A;„P,-fc(k) [{v{uj + v,k)- z/(z/, k))e + (S(cj + v,k)- S(z/, k))] 

= -u;A„jA:(0,t^; k, I/; -k,u; - z/). (94) 

Thus, a vertex correction (i.e., right-hand side of Eq fl94l) ) can arise provided that the renor- 
malized viscosity and/or the mass term are frequency dependent functions, according to the 
left-hand side of this relation. However, recall that the reference system is being subject to 
arbitrary rectilinear accelerations, and this noninertial bulk motion does introduce an ex- 
plicit time (and hence frequency) dependence into the system. Nevertheless, the correction 
to the vertex in this case is the specific one indicated in Eq. (|Ml) . We see that the vertex 
correction that enters into this relation is the one for which one of the two velocity legs car- 
ries zero wavevector but finite frequency a;, i.e., this leg corresponds to a time-dependent 
velocity zero-mode. 

Next, we come to the implication of Galilean invariance. This can be had by either 

or by taking the zero frequency limit — directly in 



inserting Eqs. (I92][93D into Eq.(|72 
Eq.l lMI) . Either way, we obtain 



t 

kn-^ [Hi', k)k^ + S(z/, k)] Pjkik) = -A„,fe(0, 0; k, u; -k, -u). (95) 

n 

As stated in [7], as a model of stationary forced turbulence, neither the viscosity nor the 
mass term will depend on frequency, so that for stationary random forcing, the constraint 
Eq. fl95l) implies that 

A„,fc(0, 0; k, z/; -k, -z/) = 0. (96) 

This specific vertex correction is in fact zero, and it corresponds to the situation in which 
one of the two velocity legs carries zero wavevector k = and zero frequency u = 0, and 
this corresponds to a velocity zero mode. 

These considerations are important. There has been a longstanding question about the 



range of validity o: 



the vertex Ward identity and the vertex non-renormalization property it 



imphes {2, s], ^, y, 3, 11, Q, Q. In jl^, we made a definitive resolution of this problem by 
discovering that the dynamic functional for Navier-Stokes theory possesses an underlying 
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and fundamental BRS symmetry. Since the fields of the BRS symmetry involve only the zero 
mode, that demonstration makes explicit that the associated vertex Ward identity is valid 
only at exactly zero momentum transfer k = 0, and not in some small limiting region (i.e., 
for — > 0) around this. The statement that vertex renormalization is not constrained by 
Galilean invariance was put forward somewhat earlier in 12l| . using physical arguments and 



employing the Reynolds decomposition. Our derivation in [13| validates the assertions made 
in [l2| through a mathematically precise and complete construction (see statements in Ref. 
jsl concerning the mathematical rigor of the path integral formalism). Except for the zero 
mode, neither Galilean nor extended Galilean invariance constrain vertex renormalization. 



VI. DISCUSSION 



It has been known for a long time that the randomly forced Navier-Stokes (NSE) equation 
can be cast in terms of a path integral 71, lot] . The underlying symmetries of 

this equation can then be treated at the level of the functional integral. In particular, the 
Galilean invariance of the NSE is formally analogous to global gauge invariance in quantum 
field theory. The well established methods for gauge fixing can then be brought to bear and 
used to restrict the sum over histories of Navier-Stokes dynamics to pick out one inertial 
frame. This gauge fixing was in fact used recently to eliminate an infinite number of spurious 
correlation functions that are implied by the non-gauge fixed functional ISj. In addition, the 
evaluation of non-Galilean invariant quantities (e.g., n-point correlations of velocity fields) 
in the standard path integral Eq.([T]) leads to over-counting of configurations and spurious 
relations, so we must fix the gauge, that is, break the Galilean invariance. If in the future, 
simulations are carried out using such functionals, our gauge-fixing procedure could possibly 
provide some numerical advantages in allowing for a faster and more efficient code. However, 



once the functional has been gauge-fixed, it is crucial to ascertain whether the wel 
a„d oft-*ed Wa.d .den«.y to. .he .ve.e ..po.e ^d ve.ex funCio. MM 



known 
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12 



remains valid or not. Regarding this question, it is important to recognize that the gauge- 
fixed theory possesses a Becchi-Rouet-Stora (BRS) symmetry [2^ which when brought out 
and exposed, restores a type of hidden Galilean invariance. The celebrated Ward identity 
is therefore re-established, but it is now understood to be a direct consequence of this BRS 
invariance, and not of the original Galilean invariance. 
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We emphasize that our BRS symmetry follows from gauge fixing and should not be 
confused with other kinds of BRS invariance, such as those treated for example, in chapter 
16 of Zinn- Justin's book 8|]. In these latter cases, it is the stochastic differential equation 
itself that is regarded as a local constraint equation for the field variable. The determinant 
of the Jacobian of this stochastic equation is then expressed as a Grassmann integral over 
ghost fields in the path integral representation of the generating functional for the correlation 
functions. An effective action is then identified, and it is this action that is shown to possess 
a simple kind of BRS symmetry, and sometimes even a kind of supersymmetry. This differs 
distinctively from our treatment here, where by contrast, the BRS symmetry emerges from 
an attempt to correct an infinity in the NSE dynamic functional arising from the underlying 
Galilean (and extended Galilean) invariance of the theory. The ghost and BRS formalism 
provides the natural language in which to establish the Ward identities. We are not required 
to "compensate" any degrees of freedom as in quantum gauge theories. Our use of the ST 
identities in the stochastic field theory of randomly stirred flows has an entirely different 
significance to the one encountered in quantum field theories. 

The NSE is invariant as well under rectilinear but otherwise arbitrary frame accelera- 
tions. This extended Galilean invariance (EGl) has been considered previously in various 



differing contexts [2^, |23|, |2J, |25|]. Just as for Galilean invariance, EGI can be regarded as 



a gauge symmetry at the level of the path integral. This higher gauge symmetry can be 
fixed, which corresponds to choosing one rectilinearly accelerating reference frame. Just as 
before, the gauge fixed theory possesses an underlying BRS symmetry which when invoked, 
leads to a generalized Ward identity for the inverse response and vertex functions. This 
new Ward identity reduces to the standard one in the limit of zero frame acceleration. It 
contains additional information regarding the vertex renormalization implied by arbitrary 
rectilinear frame accelerations, that is not provided by the "standard" Ward identity. Ex- 
;ended Galilean invariance is of interest in its own right, for by the principle of equivalence 
291], it corresponds to the invariance of the stirred flows under arbitrary time- dependent 
unidirectional background gravitational fields. 

The standard Ward identity has been appealed to on numerous occasions to make state- 
ments about the nonrenormalization of the advective or inertial term in the Navier-Stokes 
equation 
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121 ]. Recently, McComb claimed that Galilean invariance does 



not at all constrain the vertex renormalization 12(]. To quote his paper verbatim, "Galilean 
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invariance has been used as the justification of Ward identities, which in turn lead to the 
conclusion that in the perturbative renormalization group (RG) the vertex is not renormal- 
ized". Nevertheless, it is argued in Ref.[12] that "vertex renormalization is not constrained 
by this (Galilean) symmetry". Section |V] of this paper spells out explicitly what both 
the "standard" and the generalized Ward identities imply for the corrections to the bare, 
tree-level vertex. We emphasize there how our results, derived from field theory methods, 
complement, validate and confirm the assertions made in Ref [12]. To test the claim put 
forward in [l^, one must go back to the Ward identity, as well as its generalized version, 
and examine carefully what constraints are actually imposed by them. The field-theoretic 
formalism employed here is needed to work out the consequences of breaking both Galilean 
invariance and its extension. The final goal of this is to understand the vertex renormaliza- 
tion problem. The constraints on the vertex that follow from Galilean invariance and from 
extended Galilean invariance (more precisely, that follow from the respective BRS symme- 
tries of the gauge-fixed theories) are quite weak. In the case of Galilean invariance, the 
vertex correction that "couples" to the zero mode of the system is constrained to be zero. 
For EGI, the vertex correction that involves the system's bulk acceleration is computable 
in terms of the frequency dependent viscosity and mass terms of the renormalized inverse 
response function. Galilean invariance and EGI constrain only the spatial zero mode of the 
vertex but none of the higher wavenumber modes. This information is new and impacts 
directly on the physics of randomly stirred flows. These points can be appreciated from the 
diagrammatic comparison of the corresponding identities from QED and the field theory of 
the randomly stirred NSE in Sec IIVI 

It is possible to use the functional Eq([T]) without explicitly fixing the zero-mode for 
deriving perturbative expansions, as has been done up to now. Thus by implicitly fixing 
the zero mode, the perturbation expansions from the functional are consistent with those 
based directly on the NS equation itself. However, in looking for exact relations between 
correlation functions from the functional, unless one is careful, the problem of the zero 
mode then shows up and can lead to incorrect results. This mistake precisely underlies 
the confusion that has been generated over three decades in regards vertex renormalization 
claims. We have demonstrated here, following on from ISj, how to explicitly fix the zero 
mode via the analogy we identified with gauge fixing, and thus obtain a formally well defined 
functional. 
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Although we focussed attention here exclusively on the randomly forced Navier-Stokes 
equation, the gauge fixing procedure in this paper can be applied to other Galilean invariant 
theories such as the KPZ equation of random surface growth [sol, magnetohydrodynamics 



and the stochastic Burgers equation 



3l| . In 3l|, a saddle point approximation is applied 



to the path integral for the Burgers equation to calculate the tails of the probability density 
function for the velocity. A special feature of that problem is a symmetry of gauge-invariance 
type. In fact, the Burgers action is seen to be invariant under an extended Galilean trans- 
formation. This gauge de gree of freedom is then integrated over by the standard Faddeev- 
Popov trick. The aim in 311] was to show that the fluctuations around the instanton are 
free from infrared divergences. Once the action for the Burgers equation is gauge flxed and 
its EGI broken in this way, then there should be no spurious correlators there either. We 
moreover conjecture that there is a corresponding BRS symmetry that restores the EGI of 
the gauge-flxed Burgers dynamic functional with a subsequent Ward identity following as a 
consequence of this symmetry. 

A few comments regarding our use of gauge-fixing terminology are in order. With respect 
to the terminology as employed in (quantum) field theory, a local symmetry is one in which 
the associated field transformation involves parameters that depend on the spatial coordinate 
X. This is not the case for the Galilean transformation Eqs. (jSHSj), and so it is correct 
to classify this as a global (or rigid) transformation. The Galilean transformation is a 
global space and time transformation between two separate inertial frames of reference. 



In standard field theory texts (see e.g., [Sj] and [10|) the local transformations are also 



denoted as gauge transformations. Even though Galilean invariance is a global symmetry, 
we find it to useful to employ the terminology of gauge invariance and gauge fixing in this 
paper, especially since our breaking of this global symmetry can be handled with exactly 
the same methods (FP determinant, BRS, etc), as employed in gauge theories. Moreover, 
for extended Galilean invariance, the associated transformation is indeed local in time and 
so does resemble more closely a gauge theory. Note also that the convective derivative 
in the NS equation is invariant under a Galilean transformation, see e.g., the left hand 
side of Eq.f l48p . This reminds one of the invariance of the U(l) covariant derivative under 
a gauge transformation, where for EGI, the function A(t) plays the role of the arbitrary 
gauge function of QED. Since the steps for symmetry breaking in the path integral very 
closely resemble the methods employed for gauge fixing (see Appendix A), and since a BRS 
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invariance results with all its implications for Ward identities, we feel that it is useful and 
instructive to use the terminology of broken gauge invariance in this paper, albeit perhaps 
in a slightly loose manner. 



In summary this work and 13| have uncovered a new symmetry in the stirred NSE and for 
related Galilean invariant stochastic equations. Exploiting this symmetry may be useful in 
performing Monte Carlo simulations of the path integrals and actions corresponding to these 



systems 
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32|. Moreover, given the deep significance of symmetries in physical problems. 



other applications may benefit from recognizing this invariance. 
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APPENDIX A: GAUGE FIXING 

The steps outlined here closely parallel those used for gauge fixing in quantum field 



theories |26|] and are adapted to the case at hand. We first introduce the functional Aj[V] 
by the following equation: 

1 = A;[V] J dc6m% (Al) 

and denotes the result of the Galilean transformation (GT) applied to the velocity field 
V; see Eq.(|n]), and the integral is over all constant boost velocities. We assume that the 
equation f(V'^) = has exactly one solution c for any initial configuration V. If we apply 
a second GT, i.e. we next transform to a new inertial frame (double-prime) moving with a 
velocity b with respect to the primed frame, then the transformation rules Eqs.([3]|6]) tell us 
that x' = x" + ht' and t' = t" and 

(V^)^x,t) = (V'(x',^')+c)^ 

= V"(x'',0 + c + b, 

= V"(x - (c + h)t, t) + c + h = (V'')^(x, t) = V'^^. (A2) 
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This exercise is needed to prove that the Fadeev-Popov (FP) determinant A/[V] is invariant. 
The invariance of the FP determinant is needed in turn, to ensure that the volume of the 



"gauge-group" orbit can be factored out from the functional (see, e.g., [33j). The first 
requirement follows from Eq. flAip using Eq. flA2p : 

Aj^iV'^] = j d'c6%{{(V''y)) = j c/3(c + b)53^f(V^+^)) = A7^[V]. (A3) 



26|, we can prove that the volume of the gauge 



Then, repeating the standard manipulations 
group indeed factorizes out to produce an overall infinite constant. Insert Eq. flAip into the 
functional Eq.([T]) apply a Galilean transformation, then using Eq. (]A3P we obtain 

[j dc) j [W][D^]A^[V]5[f(V)]exp{-5[V,^]}. (A4) 



We still need to actually calculate the FP determinant. To calculate it, recall j33| it is 
sufficient to do so for infinitesimal "gauge" -transformations (so, infinitesimal GT's): 

A^[V] = det^^^^p|e=o = det (^)|c=o. (A5) 

We remark that the FP determinant Eq. flASp is an ordinary discrete matrix determinant: it 
is not a functional determinant. This is because the gauge parameter (= boost velocity) is 
an ordinary constant vector, not a field. Next, we consider the class of "gauge" conditions 
of the form f (V) — U(x, t) = 0, where U(x, t) is an arbitrary function. The FP determinant 
is as before because U(x, t) is unaffected by a GT. We make use of this feature to replace 
the delta function in Eq. flA4p by some other function(al) which may be more convenient for 
practical calculations. So, in this gauge the path integral Eq. (lA4p becomes 

(jdc) I [DV][DcT]Af[V]6[i{V)-Ui^,t)]exp{-S[V,a]}. (A6) 

This expression is independent of U(x, t), so we can integrate over an arbitrary weight 
functional. As usual, a popular choice is the exponential [261]: 

^[U] = exp ( - ^ y f/xrftU2(x,t)), (AT) 

for real parameter ^ > 0. For the final step, integrate Eq. (]A6P over VU to obtain 

I [DV][Da]Af[V]exp{-S[y,a] - dxdt (f (V))^}. (A8) 

This is the gauge fixed dynamic functional for randomly stirred incompressible fiuids ex- 
pressed for the gauge function f . 
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APPENDIX B: EXPANSION FOR T 



We write out the first few terms of the functional Taylor series for the effective action in 
Eq. fl7r|) . We display only those terms actually needed to derive the Ward identity in Eq. 
(125]) in this paper. In wavevector and frequency space these are: 

r[V'=',ai = ■■■ j dqiduji j rfq2t/u;2V;^'(qi,^i)fx^'(q2,cc^2)C'Vqi,-^i;-q2,-^2) 

+ dfiidui j dq2duJ2 j c/qgrfcug ^^'^'(qi, u;i)Vg^'(q2, t^2)o-^'(q3, tus) X 

ra/37 (-qi, -^i; -q2, -^2; -qs, -^3) 

+ ... (Bl) 
APPENDIX C: SPACE-TIME BOX 

The action and dynamic functional can be regularized by enclosing the system in a spa- 
tially and temporally bounded space-time box. This regularization admits a corresponding 
discrete version of Galilean invariance, implying the box Ward identity derived above in Sec 
III El from which we can demonstrate the appearance of spurious relations, in complete par- 
allel to the continuum case treated in [131]. Our box is defined from coordinates x = — L, 
y = — L, and z = — L and time defined from t = — T, subject to periodic boundary 
conditions in space and time. 

We expand the fields as 

T;c,(x,t) = ^ f j ^Va(k,u)exp[i'^{niXi+n2X2 + n3X3) -i'^jt], (CI) 

(j^ \ 3 -j^ 2y|- 27]" 

-j — dra{k,uj) exp[i— {niXi+n2X2 + 713X3) -i—jt], (C2) 



and 



(1 \ ~ 271 
-j D^p{k)exp[i—{niXi+n2X2 + n3X3)]. (C3) 
11 ^ 

The following integrals over the box are used below: 

27r 

dt exp[-z— jt] = TS,,o, (C4) 
27r 

/ rf^x exp[2— n ■ x] = L^(5n,o- (C5) 
Jo L 
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In k, uj space we will express the velocity and auxiliary fields in terms of dimensionless 
quantities as 

u„(k,t^) = t;c,(n, j)L^, (C6) 
cr„(k,t^) = cr^(n, j)— , (C7) 

and 

D,;3(k) = D„^(n)^. (C8) 

For purposes of dimensional counting, let the spatial dimensions be expressed as d{x) = L 
and temporal as d{t) = T, which also means that d{k) = L^^ and d{u!) = T~^. The 
action is dimensionless ^(5*) = 0. Then from examining individual terms in the action 
we find ^(^(k, cu)) = L^, (i(o"(k, cj)) = T^/L and d{v) = LF'/T. In configuration space 
d(w(x,t)) = L/T and d{o{^,t)) = T/L^. 

Thus the action Eq.Q in discrete wavenumber and frequency coordinates is: 



- ^X]X]cra(-n, -j)[(-27rzj + z/ra^)w„(n, j) 

- 2'KM^p^{Ti)^^vp{n-m,j -l)v^{m,l)] , (C9) 



m I 



where we defined the dimensionless viscosity i/ as z/ = uq{2tiYT / L"^ and 

Ma/37(n) = ^^^"^7(11) + n^PaM)] > (CIO) 



where 



= Saf^ - (Cll) 



In the equation Eq. ( ICQI) . everything is dimensionless. 

The Fourier transform of the Galilean transformation between prime and unprimed 
frames, analogous to Eqs.dMHD, and expressed in discrete dimensionless coordinates is 



n' 



j' = J - c- n. 



Vain,j) = v'^{n, j - c-n) + Ca5^{n)S{j), 

aa{n,j) = (T^(n,j-c-n), (C12) 
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where we defined the dimensionless boost velocity c as c = cL/T. In general c ■ n must be 
an integer for all integers n, which thus means all three coordinates of the discretized boost 
velocity c must be integers. 

With the discrete Galilian transformations Eq. flC12p in hand we demonstrate that the 
discrete action Eq. (1C9P is invariant. Begin with the term 

Gal— Trans 

m I 

^ I / . , — r n\;/,_ \ 

m 



m I 

X^X^^M^ -l-c[n- m])f;:^(m,/ - c 

I 

Therefore, we deduce the Galilean transformation of the following term in Eq. (]C9p : 



m I 

+ c^v'JnJ - c • n) + cpv'jnj - c ■ n) + CfjC^S{n)6{j). (C13) 



27ri ^ ^ o-a(-n, -j)M^p^{n) ^ ^ vp{n -m,j - l)v^{m, I) 4^''''"' 

n j ml 

= ^^^^Yl ~3 + c- n)M;^^(n) ( ^ v'p{n - m, j - / - c ■ [n - m])v'^{m, / - c ■ m) 

n j ml 

+ c^/^inj - c ■ n) + C(3v'^{n,j - c • n) + c^c^5(n)5(j)^ 

n j' m I' 

+ 27r5^5^a;(-n,-/)[c-nK(n,/). (C14) 

n j' 

In arriving at the last equality, we have used fluid incompressibility n • v = 0, the fact that 
M^^^(n)(5(n) = as well as the change of discrete summation variables j' = j — c ■ n and 
/' = / — c ■ m. Next, consider the Galilean transformation of the "propagator" term in 
Eq.dCSD: 



-i ^ ^a^(-n, - j + c • n)(-27rij + un'^)[v'^{n, j - c ■ n) + CaS{n)6{j)] 

n j 

-i ^ ^ (T^(-n, -J + c ■ n)(-27rzj + un'^)v'^{n,j - c ■ n) 

n j 

-^Xl -/)(-2vrij' + un^)v'^{n, f) - 27r^ ^a;(-n, -j')[c ■ n]v'^{n, f) . 

j' n j' 

(C15) 
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Lastly, the transformation of the "noise" term in Eq. (lC9P : 

1 / ■\ / \ / ■\ Gal-Trans 

2 Z^2^^"(-n,-j)^a/3(n)(T/3(n,j) 

n j 

= l^Yl + c ■ n)/^^^(n)(T;^(n, J - c • n) 

n j 
n j' 

Adding up the transformed terms Eqs. flC14IIC15|C16l) proves that the discrete action is 
indeed invariant under the discrete Gahlean transformation Eq (]C12p . 
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